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Abstract 

The low-field (below 5 Oe) ac and dc magnetic response of a magnetic fluid [MF] sample in the range of 305 to 
360 K and 410 to 455 K was experimentally and theoretically investigated. We found a systematic deviation of 
Curie's law, which predicts a linear temperature dependence of inverse initial susceptibility in the range of our 
investigation. This finding, as we hypothesized, is due to the onset of a second-order-like cluster-to-monomer 
transition with a critical exponent which is equal to 0.50. The susceptibility data were well fitted by a modified 
Langevin function, in which cluster dissociation into monomers, at the critical temperature [T*], was included. In 
the ac experiments, we found that T* was reducing from 381.8 to 380.4 K as the frequency of the applied field 
increases from 123 to 173 Hz. In addition, our ac experiments confirm that only monomers respond for the 
magnetic behavior of the MF sample above 7*. Furthermore, our Monte Carlo simulation and analytical results 
support the hypothesis of a thermal-assisted dissociation of chain-like structures. 
PACS: 75.75.-C; 75.30.Kz; 75.30.Cr. 
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Introduction 

The interest in magnetic fluids [MFs] has increased 
enormously in the last decade, particularly due to the 
opportunities they provide for applications in the medi- 
cal field [1-7]. Among others, MFs have been used as an 
excellent material platform for the development of mag- 
netic immunoassay [1,2], contrast agents for magnetic 
resonance imaging [3,4], and material devices for mag- 
netohyperthermia [5-7]. Deep understanding of mag- 
netic susceptibility, however, is a key issue not only 
from the fundamental point of view, but also while tai- 
loring nanosized magnetic materials for medical applica- 
tions [1-7]. The design of nanosized magnetic particles, 
taking into account the maximization of the materials' 
response in terms of their use for diagnosis, imaging, 
and therapy, requires the knowledge of the temperature 
dependence of the magnetic susceptibility under the 
action of applied dc and ac fields. 
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In this context, the widely accepted concept is a linear 
relationship between the inverse initial magnetic sus- 
ceptibility (11%) and the temperature {T), which is 
accounted for by the first-order Langevin function. 
Nevertheless, unusual deviations of linearity at tempera- 
tures within the range of interest for the medical appli- 
cations, with no conclusive explanation yet, have been 
reported [8]. Since an interaction among particles in a 
MF sample, either modulated or not by external fields, 
cannot be ignored and leads to expontaneous agglom- 
eration in clusters or chain-like structures (dimers, tri- 
mers, etc.) [9-11], the superlinear deviation of the \l% 
versus T curve has been attributed to magnetic dipolar 
interaction among the nanosized particles [12]. The 
Langevin function, however, includes no interaction 
among the suspended particles, and therefore, clusters 
are ruled out from the classical description. Neverthe- 
less, at high enough temperatures, the phenomenon of 
thermally assisted cluster disruption within MF samples 
has been reported [13]. Therefore, to understand the 
underlying physics of this nonlinearity in magnetization, 
a more complete physical model is highly demanded. 
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The new model should take into account both mono- 
mers and clusters in MFs. 

In this study, we report the unusual superlinear devia- 
tion in the temperature dependence of the inverse initial 
magnetic susceptibility in a magnetite-based (Fe 3 0 4 ) MF 
sample in the temperature range of 305 to 360 K. The 
observed breakdown of Curie's law, which scales linearly 
the inverse susceptibility with temperature, indicates 
that besides the usual tendency of alignment of mag- 
netic moments with the applied field, there exists an 
additional thermally assisted physical process connected 
to cluster disruption within MF samples. Therefore, we 
propose a model in which a chain-like disruption at a 
typical transition temperature (T*) is incorporated. We 
found that the extended model reproduces quite well 
the observed superlinear deviation. Additionally, the 
above-mentioned superlinear deviation of the (1/^) x T 
data for chains of particles is supported by Monte Carlo 
[MC] simulation and herein incorporated. 

Experimental description 

In order to verify the universality of the experiment 
regarding the relationship of inverse initial susceptibility 
and temperature, we explored the experiments under dc 
and ac magnetic fields, the latter at different frequencies. 
The MF sample used in our experiment was a commer- 
cial magnetic colloid (EFH1, Ferrotec Corporation, 
Santa Clara, CA, USA), consisting of Fe 3 0 4 magnetite 
nanoparticles (mean particle diameter of 10 nm) sus- 
pended in light mineral oil. The applied magnetic fields 
were 5 and 2 Oe (amplitude) for dc and ac experiments, 
respectively. The ac experiments were performed at 123 
and 173 Hz. 

The chord magnetic susceptibility (% = MIH) can be 
obtained by measuring the sample's magnetization (M) 
induced by a weak external dc/ac field (H). The home- 
made experimental system designed to measure the MF 
sample's initial susceptibility is schematically shown in 
Figure 1. The Helmholtz coils were driven by a dc/ac 
power supply to generate a uniform dc/ac excitation 
field at the sample's position; the direction of the dc 
applied field can be altered by a relay. The temperature 
of the MF sample was measured using a thermocouple 
(PtlOO) and converted using a temperature transmitter 
(RS1852290). A weak magnetization signal collected 
from the MF sample was detected using the two-axis 
giant magnetoresistive sensors (HMC1022) with an 
amplification circuit. Both temperature and magnetiza- 
tion data of the sample were acquired through a data 
acquisition card (PCI-6251, National Instruments, Aus- 
tin, TX, USA). Lab VIEW was used as a human-compu- 
ter interface for observation and control. 

The experiment includes three steps: Firstly, the sam- 
ple is placed inside a water bath, slowly heated up to 



the water's boiling point temperature, and kept there for 
a certain period of time. While keeping the sample 
holder stably hot, a weak uniform dc/ac field was 
applied to the MF sample. The second step is to record 
the inverse susceptibility (x 1 = HIM) versus tempera- 
ture curves. The heating circuit is shut off, allowing the 
sample of cooling down naturally while synchronously 
recording both the temperature and the chord suscept- 
ibility. To avoid stray magnetic field interference, when 
the dc field is applied, the direction was set to invert 
every 0.5 s. The magnetization was evaluated according 
to the difference between the two measurements, before 
and after switching the relay. For ac susceptibility mea- 
surement, the digital phase-sensitive detection algorithm 
was used to calculate the amplitude of the magnetiza- 
tion. The last step is data processing. Regarding the data 
obtained from repeated experiments at different excita- 
tion fields, we took a temperature point every 0.5 K in 
the temperature range of the experiments, and the mea- 
sured values at temperature points within a permissible 
measurement error range are considered to be achieved 
at the same temperature. Accordingly, the mean values 
of repeated experiments were obtained to draw the 
curves of inverse susceptibility versus temperature. Typi- 
cal {11 X) x T data are shown in open symbols in Figures 
2 and 3. Note the superlinear behaviors of the (1/x) x T 
data at the higher temperature end in both Figures 2 
and 3. 

Model and discussion 

The susceptibility model for pure monomers 

Considering only isolated particles (monomers) with no 
interparticle interaction, the MF sample's magnetization 
is well described by the Langevin function: 

M = <pM s (coth{mH/kT)-kT/mH), (1) 

where cf) is the concentration of the nanoparticle (the 
number of particles per unit volume), m is the nanopar- 
ticle's magnetic moment, H is the applied field, k is the 
Boltzmann constant, and T is the absolute temperature. 
The nanoparticle's magnetic moment (m) is described in 
terms of the saturation magnetization (M s ) as m = M S V, 
where V is the nanoparticle's volume. At weak enough 
applied fields, the sample's susceptibility^ = dM/dH 
approaches the chord susceptibility % = MIH. Then, the 
inverse initial susceptibility can be written as: 

X' 1 = (H/<p) /[M s (coth(mH/fcT) - kT/mH)]. (2) 

Equation 2 describes the {11 X) x T data based on the 
first-order Langevin function. Note that the dashed line 
in Figure 2 is the curve fitting of the data (open sym- 
bols) using Equation 2, showing the expected linearity 
in the temperature range of interest. However, our 
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Figure 1 Schematic diagram of the measurement system suitable for dc/ac initial magnetic susceptibility detection. 



experimental data (open symbols in both Figures 2 and 
3) revealed a superlinear trend at the high temperature 
end, making the fitting procedure using Equation 2 visi- 
bly poor (see inset (a) of Figure 2). We then hypothe- 
sized that the assumption of only monomers in the MF 
sample in the temperature range of 305 to 360 K no 
longer holds. 

Discussion of the susceptibility model including 
dimers 

The explanation of the superlinear behavior observed on 
the (1/%) x T experimental data displayed in Figures 2 
and 3 (open symbols), as we claim, is due to the pre- 
sence of a fraction of chain-like structures within the 
MF sample, in addition to monomers, more likely 
dimers [13,14], and the temperature dependence of the 



dimer fraction in a critical way around T*, as discussed 
below. 

The probability of agglomeration and disruption of 
suspended nanoparticles in MFs is assumed to be 
dependent upon the relative strength of magnetic, van 
de Waals, electrostatic, and steric interactions and ther- 
mal energy. The literature [15-17] describes that mag- 
netic dipolar interaction held particles together, whereas 
electrostatic interaction and thermal energy work 
together taking nearby particles apart. It can be inferred 
from the experimental observations that there is a 
dynamic balance between the relative content of mono- 
mers and dimers within the simplest model picture of a 
magnetically textured MF sample. This balance will be 
broken when there is a change on the parameters gov- 
erning the energy terms involved, and the relative 
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Figure 2 Experimental data and the fitting curves using the Langevin function. Experimental data under a 5-Oe dc external field (open 
squares) and the fitting curves using the first-order Langevin function (dashed line) and the modified Langevin function (solid line). Inset (a) is a 
zooming in the temperature range of 350 to 360 K. Inset (b) is the temperature dependence of P 2 . 
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content of monomers and dimers will change accord- 
ingly. As the thermal energy increases, the dimers tend 
to disrupt into monomers, leading to a decrease of the 
dimers' relative content. Inversely, within this simplest 
model picture, when lowering the temperature of a MF 
sample, the monomers tend to agglomerate into dimers, 
leading to a decrease of the monomers' relative content. 

The already described process of the thermal-assisted 
dimer disruption within a MF sample [12] is related to a 
second-order phase transition at a critical temperature 
T*. Below T*, the suspended nanoparticles are found as 
monomers and dimers, whereas above T*, nanoparticles 
are essentially isolated within the MF sample while the 
relative content of dimers drops down critically to 
nearly zero around T*. Below the critical temperature 
T*, monomers and dimers do coexist within the MF 
sample. In order to describe this mixed (magnetically 
textured) system and the influence of the temperature 
upon it, we have modified Equation 1 including two 
contributions, namely terms due to monomers and 
dimers. Below T*, the probability of finding monomers 
(Pj) and dimers (P 2 ) within the MF sample scales with 
the temperature according to: 

Pi = T/r &P 2 = (l - T/T*y. (3) 

Considering the classical Landau's theory of second- 
order phase transition, f} in Equation 3 represents the 
critical exponent (/3 = 0.5) as long as P 2 properly 
describes the order parameter of the thermally assisted 
dimer disruption process. Thus, the modified Langevin 
function describing the actual magnetization of the MF 
sample is: 

M - ipP,M sl [coth^nH/kr) - [Kr/m 1 H)]-HpP 2 M S 2 [coth(m 2 H/hT) - (feT/m 2 H)] , (4) 

where m.\ and m 2 are the magnetic moments of 
monomers and dimers, respectively. In Equation 4, (\>P\ 
and 4>P 2 represent the content of monomers and dimers, 
respectively. For abbreviation, we write Ci = m\HlkT 
and £2 = miHIkT. According to Equations 3 and 4, the 
inverse susceptibility describing the magnetically tex- 
tured MF sample can be written as: 

PT 1 = H/j^dMs, [coth(£,) - (1/50] + <pP 2 M s2 [cothfe) - (1/fe)]). (5) 

The low-field dc susceptibility data are shown in Fig- 
ure 2. The solid line represents (see Figure 2) the best 
curve fitting of the (llx) x T experimental data (open 
squares) using Equation 5. Included here for compari- 
son, the dashed line in Figure 2 represents the fitting of 
the data using Equation 2. Note that ferromagnetic reso- 
nance data (resonance line splitting) have been used to 
describe dimer disruption in a nickel ferrite-based ionic 
MF sample, providing values of T* (340 K) and ft (0.42), 



the latter in reasonable agreement with the value we 
found in the present study [13]. Our data indicate that 
above T*, the suspended particles within the MF sample 
investigated are essentially isolated (monomers). The 
inset (a) of Figure 2 shows a detail of the (1/^) x T data, 
emphasizing the superlinear behavior at temperatures 
close to the typical dimer disruption temperature (T*). 
The inset (b) of Figure 2 shows the temperature depen- 
dence of P 2 , here, describing the order parameter asso- 
ciated to a second-order-like phase transition. 

Experiments on low-field ac susceptibility at different 
frequencies (123 and 173 Hz) are shown (symbols) in 
Figure 3. The saturation magnetization of monomers 
(M S i) and dimers (-M S2 ) in the same MF sample was 
considered to be constant. Af si and M S2 were firstly 
acquired by the fitting of any applied field (dc or ac) 
and frequency, and then both iVf sl and M S2 were used as 
known parameters to fit different data sets. We found 
different susceptibility responses while changing the fre- 
quency of the applied field. The fittings of the experi- 
mental data (solid lines) using Equation 5 show that the 
critical temperature decreases from T s = 381.8 down to 
T" = 380.4 as the frequency of the ac excitation field 
increases from 123 to 173 Hz, in agreement with 
recently reported results [18]. 

Monte Carlo simulation of second-order-like dimer- 
monomer transition 

Because a dipolar interaction between neighboring mag- 
netic nanoparticles favors parallel alignment of their 
magnetic moments, a strong interaction may induce 
coalescence of nanoparticles and formation of magneti- 
cally coherent clusters such as dimers. Then, monomers 
and dimers may coexist in the system with an inter- 
mediate concentration of nanoparticles. This model pic- 
ture allows us to describe the system using the Ising 
model in which we consider AL identical magnetic 
nanoparticles disposed in one-dimensional chains. We 
assume that each nanoparticle can be represented by a 
magnetic monodomain, with a magnetic moment M s cr, 
along the direction of an external magnetic field, where 
Oi = + 1, i = 1, 2,..., N p . Then, the total energy of the 
nanoparticle assembly is given by 

E = -J 52 . \ Oiaui - h Y\ . *. Oi, (6) 
^— 't=i *—'i=\ 

where the sum is extended to the nearest neighbors 
only, / is the nearest neighbor exchange integral, and h= 
M S H. The equilibrium magnetic moment configuration 
is obtained by a MC simulation using the standard 
Metropolis algorithm. It is interesting to find out that 
for systems such as MFs in which suspended nanoparti- 
cles are found as monomers and dimers, the 
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magnetization and susceptibility can be calculated analy- 
tically through the partition function of the system, 
Z = J] g- fi n/ feT , where E n is the energy of the system in 
the «th configuration. For instance, we can use the fol- 
lowing four configurations in order to describe dimers | 
CTitT 2 >=|++>.|+->,|-+>,|->- 

Temperature-dependent disruption process of dimers, 
i.e., variation on the number of dimers (AT d ) in MFs, can 
be modeled by a parameter y, which is defined as 
Y = N d /(N m + N d ) = ^l-(T-Ti)/(T f -Ti)/2, where 
N m is the number of monomers, AL = N m +2N d , and T l 
(X f ) is the initial (final) temperature of the simulation. 
Notice that y describes the system in a wide range of 
configurations in regard to the disruption process, for 
instance, from a MF with the same number of mono- 
mers and dimers (N d = N m ) at T = T t to a configuration 
in which dimers are completely disrupted (N d = 0) at T 
= Tf. Then, the system's susceptibility consisting of 
monomers and dimers (complex system) is given by 

X = N p [(1 - Y ) Xm + YXdi/a + Y). (7) 

where % m {%&) is the susceptibility of the monomers 
(dimers). Figure 3 shows the temperature dependence of 
the inverse susceptibility obtained by MC simulation for 
monomers (dashed line), dimers (dash-dotted line), and 
their complex (solid line) at \iHIJ = 0.1, kT t IJ = 0.5, and 



kTjIJ = 10. Notice that the inverse susceptibility of 
monomers is larger than that of dimers, and both of 
them show almost linear dependence on temperature. 
However, an unusual superlinear deviation is observed 
in the monomer-dimer complex curve. From Equation 
7, one sees that % mainly depends upon % m , % d , and y. 
Therefore, the variation of y with temperature, i.e., the 
change of the number of dimers in a MF sample, results 
in the deviation of the susceptibility of MFs from linear- 
ity. Finally, the analytical solution of the monomer- 
dimer complex system, as shown in inset of Figure 4, 
gives a further strong support for this model picture. 

The pure monomer scenario above T* 

Experimental investigation of the inverse susceptibility 
versus temperature was also performed in the tempera- 
ture range above T*. In order to prevent volatilization of 
the MF solvent (mineral oil), the highest temperature 
achieved in our experiment was restricted to 455 K. 
According to the analysis presented above, all dimers 
are assumed to be disrupted into monomers when the 
temperature of the MF sample is set above the critical 
temperature T*. This means that above T*, the relative 
content of dimers tends to be 0% whereas the relative 
content of monomers tends to be 100%. The inverse 
susceptibility-temperature curve in the temperature 
range of 410 to 455 K is plotted in Figure 5. The 
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Figure 4 Inverse susceptibility of monomers, dimers, and their complex. Inverse susceptibility of monomers (dashed lines), dimers (dot- 
dashed lines), and their complex (solid lines) as a function of temperature for fjH/J = 0.1, obtained by MC simulation. Inset presents the 
correspondent analytical solution. 
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experimental data shown in Figure 5 were fitted by the 
monomer model, herein described by: 

X -1 = H/[(pM s (coth(mH/kT) - kT/mH) + A], (8) 

Note that A in Equation 8 means just a scale factor 
for the experimental data. The solid line in Figure 5 has 
an excellent agreement with the experimental data and 
presents the best curve fitting achieved using Equation 
8. Furthermore, the monomer's saturation magnetization 
(Ms) obtained from the data recorded in the tempera- 
ture range above T* (from 410 to 455 K) is very much 
close to the value found from the data recorded in the 
temperature range below T* (from 305 to 360 K). 
Although fitting of the data presented in Figure 5 (open 
symbols) can be performed using the dimer-monomer 
model, represented by Equation 5, there is no agreement 
between the /Vf sl obtained from below and above T*. 
This indicates that dimers can be assumed to be totally 
disrupted into monomers when the temperature is 
increased above T*, and the critical temperature T* can 
be actually used to describe a MF sample. 

Conclusion 

In conclusion, the usual linear temperature dependence 
of the inverse initial susceptibility (dc or ac) of MFs at 
lower temperatures is found whereas at higher tempera- 
tures, an upward deviation of linearity is observed. This 



superlinear behavior is attributed to a thermal-assisted 
disruption of dimers into monomers which is described 
by the classical Landau's approach for second-order 
phase transitions. The experimental observations in the 
temperature ranges below and above the critical tem- 
perature T* are well fitted by the model picture pro- 
posed here, in which an extended Langevin function 
including the thermal criticality of the MF system is 
adopted below T*. Our findings are strongly supported 
by both MC simulation and analytic analysis. 
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